Input-output relations at dispersing and absorbing planar multilayers for the 
quantized electromagnetic field containing evanescent components 
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By using the Green- function concept of quantization of the electromagnetic field in dispersing 
and absorbing media, the quantized field in the presence of a dispersing and absorbing dielectric 
multilayer plate is studied. Three-dimensional input-output relations are derived for both ampli- 
tude operators in the k-space and the field operators in the coordinate space. The conditions are 
discussed, under which the input-output relations can be expressed in terms of bosonic operators. 
The theory applies to both (effectively) free fields and fields, created by active atomic sources inside 
and/or outside the plate, including also evanescent-field components. 

PACS numbers: 42.50.-p, 42.50.Ct, 42.50.-p, 42.79.-c 



I. INTRODUCTION 

Research on quantum statistical properties of light has 
been a subject of intense investigations and discussions. 
In view of the very wide-spread potential applications, 
the influence of material bodies on the quantum features 
of light must be thoroughly investigated. Typical quan- 
tum effects that arc closely related to the change of the 
quantum vacuum due to the presence of material bod- 
ies are the Casimir effect [1] (for a review, see [2]) and 
the Purcell effect [3] (see also [4]). Though there has 
been a large body of work on the problem of quantiza- 
tion of the electromagnetic field in dielectric media (for 
a review, see [5]), most work has been concentrated on 
nonabsorbing media. Roughly speaking, there have been 
two routes to treat radiation in absorbing media. In the 
first, attention has been restricted to equilibrium field 
correlation functions that are calculated by employing 
the dissipation-fluctuation theorem [6]. In the second, 
explicit field quantization has been performed (for a re- 
view, see [7]). 

A consistent quantization scheme for the electromag- 
netic field in bulk material has been given by Huttner and 
Barnett [8] . Using the mesoscopic Hopfield model [9] for 
a dielectric, they have diagonalized the bilinear Hamil- 
tonian of the system that consists of the electromagnetic 
field, a harmonic-oscillator polarization field, and an infi- 
nite set of harmonic-oscillator reservoir fields responsible 
for absorption. By starting directly with the phenomcno- 
logical Maxwell equations for the macroscopic electro- 
magnetic field, the scheme can naturally be extended 
to arbitrary inhomogeneous media [10] characterized by 
a space- and frequency-dependent complex permittivity 
that satisfies the Kramers-Kronig relations. The theory 
is based on a source-quantity representation of the elec- 
tromagnetic field, where the field operators are expressed 
in terms of a continuous set of fundamental Bosc fields 
via the Green tensor of the classical problem (for details, 
see also [7]). 

The Green-tensor formalism is well suited to study 
the behavior of the quantized electromagnetic field in 



the presence of dispersing and absorbing bodies. In par- 
ticular, it has successfully been applied to the study of 
input-output relations [11, 12], the spontaneous emission 
[13-15], the resonant energy exchange [16], the resonant 
dipole-dipole interaction [17], and the Casimir effect [18], 
and also various geometries (including multilayer struc- 
tures) have been considered. It is worth noting that the 
theory not only takes into account material absorption 
in a consistent way but it also includes automatically 
evanescent- field contributions generated, e.g., by radiat- 
ing atoms that are very close to the bodies under consid- 
eration. 

For the particular case of a dispersing and absorb- 
ing dielectric structure in vacuum, without active light 
sources at any finite distance from the structure (and 
without active light sources inside the structure), a 
ratherinvolved hybrid quantization scheme has been pro- 
posed [19] and used to study the problem of three- 
dimensional input-output relations [20, 21]. In the 
scheme, the whole space is divided into two regions, 
namely a dielectric scattering region and a vacuum re- 
gion that surrounds the scattering region. The quan- 
tization in the vacuum region is performed on the ba- 
sis of a mode decomposition of the electromagnetic-field 
operators in a similar way as in vacuum quantum elec- 
trodynamics, whereas the contribution of the dielectric 
scattering region is taken into account by applying the 
Green tensor formalism [10]. The argument is [21] that 
the Green tensor formalism alone is not complete, be- 
cause it does not explicitly treat the boundary of the 
medium with free space and thus does not describe ex- 
plicitly the scattering and emission processes which are 
usually experimentally investigated. This is of course 
not the case as the above mentioned applications [11- 
16, 18] and the following study of the three-dimensional 
input-output relations at multilayer plates clearly show. 
Moreover, when there are active light sources at finite 
distances from the scattering region, then the incoming 
field contains both propagating and evanescent compo- 
nents. The latter ones, which are typically observed for 
small distances, are not included in the hybrid quantiza- 
tion scheme [19] and the input-output relations derived 
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from it [21]. However, they are automatically included 
in the Green tensor formalism [10]. 

In the following we consider the problem of the 
three-dimensional input-output relations for the electro- 
magnetic-field operators at the boundaries of a dielectric 
planar multilayer structure in more detail, by applying 
the Green tensor formalism and extending our previous 
one-dimensional analysis [11] to three dimensions. To be 
quite general, we allow (i) for the presence of active light 
sources at arbitrary positions inside and/or outside the 
multilayer plate and thus (ii) for both propagating-field 
components and evanescent-field components. In partic- 
ular, in cavity QED the active sources are typically inside 
the plate, and the outgoing fields are not only determined 
with the incoming fields and the noise fields (unavoid- 
ably associated with material absorption) but also with 
the fields generated by the sources inside the plate. On 
the contrary, in optical near-field microscopy the active 
(probe) sources are typically outside the plate, but near 
its surface. 

Employing the well-known three-dimensional Green 
tensor for a multilayer plate [22], we introduce ampli- 
tude operators for the input and output fields as well as 
for the fields inside the plate and derive input-output re- 
lations both in the two-dimensional Fourier space and in 
the coordinate space. Finally, the problem of introduc- 
tion of bosonic input and output operators in the two- 
dimensional Fourier space is considered in detail. 

The paper is organized as follows. In Sec. II the quanti- 
zation scheme is briefly summarized and the Green tensor 
for a multilayer plate is introduced. In Sec. Ill the gen- 
erally valid input-output relations are derived, and the 
problem of formulating input-output relations for bosonic 



field operators is studied. A summary and some conclud- 
ing remarks are given in Sec. IV followed by an appendix, 
in which relevant commutation relations arc derived. 



II. BASIC EQUATIONS 

To study optical fields interacting with active sources 
in the presence of dispersing and absorbing (linear) di- 
electric bodies, we first note that on a length scale that is 
large compared with interatomic distances in the bodies, 
the effect of the bodies can be described within the frame 
of macroscopic Maxwell equations in terms of a spatially 
varying permittivity which is a complex function of fre- 
quency. This concept, which is widely used in classical 
optics also applies in quantum optics. 



A. Field quantization in media 

Let us assume that the active light sources are neu- 
tral atoms and consider an arbitrarily inhomogeneous 
medium characterized by a permittivity 

e(T,u)=e'(T,u)+ie"(T,w), (1) 
where the real part e'(r,u>) and the imaginary part 
e"(r, uj) are necessarily related to each other through the 
Kramers-Kronig relations, due to the causality princi- 
ple. The motion of the atoms and the medium-assisted 
electromagnetic field is then governed by the multipolar- 
coupling Hamiltonian [7] 
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where A numbers the atoms, and a A numbers the 
charged particles inside the A-t\i atom. Further, r aA 
and p aA are respectively the operators of coordinates 
and canonical momenta of the particles, and 

P^( r ) =y2la A (r aA -r A ) / d\5[r-r A -X(r aA -r A )] 

(3) 

is the operator of the polarization of the A-th atom at 
position r A . The f(r, uj) [and F(r, w)] are bosonic field 
operators that play the role of fundamental variables of 
the electromagnetic field and the medium, including a 
reservoir necessarily associated with material absorption. 



The commutation relations for the operators f(r, u) and 
F(r, u>) are 

[^(r, W ),/t,(r',o/)] = V^-^ (3) (r-r'), (4) 



[^(r, W ),jV(r>')] =0, 



(5) 



where the Greek letters label the Cartesian coordinates 
x, y, z. In Eq. (2), the operators B(r) and E(r) of the 
medium-assisted electromagnetic field are expressed in 
terms of the fundamental variables as follows: 



B(r) = f 

Jo 



duj B(r, u) + H.c., 



(6) 
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E(r) = / dwE(r,w) +H.c, 



B(r,w) = (iw) _1 V x E(r,w), 



E(r,w) = ifi \ — w 



/" d 3 rVe"(r»G(r, r', w) • f (r', 



(7) 



(8) 



(9) 



. active sources .... 



where the integration should be performed over all space, 
and G(r, r', oj) is the classical Green tensor, which can be 
found from the equation 

UJ 2 

V x V x G(r, r', w) - — e(r, w)G(r, r', w) = (5 (3) (r - r') 

(10) 

together with appropriate boundary conditions at infin- 
ity. Equations (6 - 9) can be considered as generalization 
of the familiar mode decomposition that would apply if 
dispersion and absorption could be disregarded. Instead 
of dealing with equations of motion for mode operators, 
equations of motion for the fields f (r, oj) must be treated. 

It should be pointed out that since the real part e'(r, us) 
and the imaginary part e"(r, w) of the permittivity are 
related to each other through the Kramers-Kronig rela- 
tions, e"(r,oj) cannot vanish identically for really exist- 
ing media. Clearly, e"(r',oj) can be very small, so that 
a/ e"(r', u>) G(r, r', u>) in Eq. (9) is very small in certain 
areas of space (r'). However, this statement says nothing 
about the magnitude of the total integral over the coor- 
dinate r', because of the following integral relation for 
the classical Green tensor 

J dV^e"(r', U )G^ (r, r', w)G£, v (r", r', w ) 

= Im [G^(r,r'»], (11) 

where we have adopted the convention of summation over 
repeated vector-component indices. In particular, this re- 
lation enables one to include also vacuum-like areas in the 
consideration. Thus, all the calculations are to be per- 
formed by assuming a permittivity close to unity with a 
small but finite imaginary part in those areas, and at the 
end the permittivity may be set equal to unity. In prac- 
tice, experimental realization of (macroscopic) vacuum 
areas is of course fictional. 



B. Multilayered planar structures 

As mentioned, the quantization scheme is valid for an 
arbitrary space dependence of the permittivity. Here, 
we consider a multilayered planar structure (Fig. 1), 
whose permittivity is defined in a stepwise fashion (the 
^-direction is chosen to be perpendicular to the layers): 



(12) 




z = 0' 



E. 



z = H 



ACn) 



O) 



FIG. 1: Scheme of the multilayer dielectric plate. The hatched 
regions indicate the presence of active light sources. 



where 



X j (z) = 



1, if z 6 jth layer, 
0, otherwise, 



(13) 



and Ej(uj) is the (complex) permittivity of the j-th layer. 
In the above, the index j labels the region on the left of 
the plate (j = 0), the region on the right of the plate 
(j = n )i an d the layers of the plate (J = 1, . . . , n — 1). 
For simplicity, we express the z-coordinate dependence 
in shifted coordinate systems, introduced in each layer 
separately, so that the range of the z-coordinate is taken 
to be — oo < z < for the region on the left of the plate 
(j = 0) , < z < oo for the region on the right of the plate 
(i = n )> an d 0< z <dj for the j-th layer of the plate with 
thickness dj (j = 1, . . . , n — 1). 

Exploiting the translational symmetry in the (a;y)-pla- 
ne, we may represent the Green tensor as a two-dimen- 
sional Fourier integral 

gW>{t,j?,u) = -^2 J&ke^i-rt G {ji '\z,z'Xu)), 

(14) 

where p = (x, y), and k = (k x , k y ) is the wave vector par- 
allel to the layers. The notations '(v, r', ui) and 
G^i ' (z, z', k, w) indicate that z varies in the j-th layer 
and z' in the j'-th layer. Inserting Eq. (14) into Eq. (9), 
we may write the electric-field operator E^ (r, lj) as a 
two-fold Fourier transform, 

E W (r,^) - J d 2 fc e ' k "E(%,k, W ), (15) 

where 

EW(z,k,w) = iujuo 
n /> 

xV / dz'G to,) (z,z',k,cj)-j^")(z',k,w)(16) 

r=0 J m 
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([j'} indicates integration over the j'-th region). Here, 
the Fourier transformed current operators 



P> (z, k, W )= W ,p e» f 0) (z, k, W ) (17) 



have been introduced, where the operators {^(z,k, u>) 
are the Fourier transforms of the original bosonic field 
operators f ^ (r, uj), 

f U) (r, lo) = J d 2 k e ik " f W) k, W ). (18) 

The Green tensor G [jj '\z, z', k, uj) taken from Ref. [22] 
may be written as 



G fo,) (z,z',k,w) 



S(z-z')+g^ 3 '\z,z',k,uj), (19) 



and tj/ji = (Pj/Pji)tj>/j and Tj/y are, respectively, the 
transmission and reflection coefficients between the layers 
j' and j. Finally, the unit vectors e^2(k) in Eqs. (21) and 
(22) are the polarization unit vectors for TE (q = s) and 
TM (q — p) waves, 



Af) 



3 p± 



(k) 



(k) 



k 

k 



1 



(27) 
(28) 



The 'propagation constant' {3j determines the propaga- 
tion behavior in z-direction of the waves in the j-th re- 
gion. Note that in case of vacuum the waves are propa- 
gating only for (ij — (3j (i.e., uj/c> k). They are evanes- 
cent for (3j =i(3j (i.e., uo/c<k). 



where 

g("\z,z',-k,w) 

q=p,s 

+ £{< (z, k, w) Z(i£(>(z', k, u>)Q{j' - j)] (20) 

(cr p = 1, a s = —1). Note that for j =j' one should write 
Q(z-z') instead of Q(j — f) [and Q(z' — z) instead of 
©(.?' - j)]- In Eq. (20), the functions (k,uj, z) and 
£ 3 ^~ (k, uj, z) denote waves of unit strength traveling, re- 
spectively, rightward and leftward in the j-th layer and 
being reflected at the boundary, 

£V»{z,-k,w) 

= e^(k)e^^- d ^ +r 9 j/n eVl(k)e- i M*- d i\ (21) 

e?<(z,-k,w) = e^(k)e-^-+r| /0 e^(k)e^-, (22) 
and 

1 tl /o e^ t^eW' 



q D qj D, 



where 



D vi = l-r q /n r q , e 2i ^ 

13 j/0 j/n 



(23) 



(24) 



(do = d n = 0). Here, 

= y/k] - A 2 =(3> j+ tf> (p>. , 0> > 0) (25) 
(k= |k|), where 



k 3 = yfe^jj) - c = k ' 3 + ik " . k'j > 0) , (26) 



III. INPUT-OUTPUT RELATIONS 

The equations of motion which follow from the Hamil- 
tonian (2) in Sec. II A together with the Green tensor 
as given in Sec. II B can be used to study various effects 
of the atom-light interaction in the presence of dispers- 
ing and absorbing multilayered planar structures (see the 
examples mentioned in Sec. I). Here, we restrict our at- 
tention to the general relation between the fields outside 
and inside a multilayer plate regardless of the generation 
of the fields by the active light sources. The theory thus 
includes the general case, where sources are present in 
the both regions, outside and inside the plate. 

If the incoming fields, incident on the two boundary 
planes of the plate (i.e., z = 0~ for j = and z = + for 
j =n; cf. Fig. 1), are known as well as the fields gener- 
ated inside the plate, one can calculate the fields outgoing 
from the two boundary planes by means of input-output 
relations. Note that any two planes z — z^ < 0~ and 
z = z( n ) > + for j = and j = n, respectively, also can 
be used in principle. As we shall see, these input-output 
relations are valid for both passive and active devices, for 
arbitrary layer materials and arbitrary media surround- 
ing the plate and for both propagating-field components 
and evanescent-field components. 



A. Input-output relations in the k-space 

Here and in the following we restrict our attention to 
the electric field, noting that the corresponding expres- 
sions for the magnetic field can readily be obtained by 
applying Eq. (8). Following the line suggested in Ref. [11] 
and substituting the Green tensor (19) for j = and j = n 
into Eq. (16), we may decompose the field operators 



■5 



E(°)(^,k,w) and E(")(z,k,w) in the form of 

E<°>(*,k, W )= £ [e^(k)EW(,,k, W ) 



and 



+ e^(k)E^ ut (z,k,^ , (29) 



E<»>(*,k, W ) = ]T [e^(k)E^(z,k,^) 



q=p,s 



+ e$(k)E<l t (z,k,a,) . (30) 



Here, the operators 



f dz'e-^'j(°)(z',k^).e( 4 ) (k), (31) 

J —OO 



P QO 

/ dz 1 e l ^ z ' f n \z' ,\h,uj) •e^(k) 



E(° ) ut (z,k^)^e^^E(° ) ut (k,c,) 



/•0 

+ e - l/3oZ / dK / e <A > z 'j( >(« / ,k,w)-e5^(k), 

J Z 



(z,k, W )=e^*EW (k, W ) 



(33) 



+ e l/3 " z / dz' e^" 2 ' jW (z', k, w) • (k), (34) 
Jo 



respectively, play the role of input and output amplitude 
(32) operators, where the input-output relations 



E 

Q ULl 

E^lt(k,^ 



^S(k >W ) 

*0/n( k ' W ) </ ( k ^)y 



r 0/n( k ><") C/ ( k ' W ) 



n-l 

£ 



^(k.w) 0?2_(k, W ) 
^ + (k, W ) <t>f n -{^) 



Eg(k,o;) 
E«(k,u,) 



!, (35) 



r 



are valid (for the commutation relations, see the ap- 
pendix). They relate the output amplitude operators at 
the boundary planes of the multilayer plate to the input 
amplitude operators at the boundary planes, 



E 



(o) 



t( k , w ) 



E 



(0) 



in, out V"-? ^ J in, out ^0 

E^, out (k,c) = Eg iOUt (^,k,a;) 
and the intraplatc amplitude operators 
Eg(k,o;) = -"" 



2 = 0" 



2 = + 



(36) 
(37) 



2(3, 



f ' dz'e^ z 'p\z',-k,w) -eg(k) 
Jo 



(38) 



(j = 1, 2, . . . , n — 1), which are associated with the ex- 
citation inside the layers of the plate. In Eq. (35), the 
(^-coefficients read 



A3) 

P q0+ 



3/0 



D 



-r q 

'j/n> 



Aj) 
D q0- 



p, 09) 

U 13 



t q , e if} i d i 
Aj) _ 3_[n 



D 



<I3 



Aj) 3/n 
J qn — 



3 to- 



rn 



It should be pointed out that the first term on the right- 
hand side in Eq. (19), which gives rise to a local contri- 
bution to the electric field, has been omitted in Eqs. (29) 
and (30). Though this contribution is irrelevant for the 
incoming and outgoing fields, it must be included in the 
overall field operator in general, even if there are effec- 
tively no sources at the points of observations. 

It is not difficult to prove that (similar to the onc- 
dimensional case [11]) the z-dependent amplitude opera- 
tors (31) - (34) obey quantum Langevin equations, 



-^j(°>(z,k,c) 
2po 



•eg(k), (41) 



+ ^j(°)(z,k,c).4°) ( k), (42) 

and similar equations are valid for E^(z,k, ui) and 
E^ ut (z,k, uj). These equations together with Eq. (35) 
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render it possible to easily calculate the input and out- 
put fields at any position outside the plate. Needless 
to say that other than the boundary planes z = CP and 
z = + of the plate can be chosen as reference planes for 
formulating the input-output relations. 

Equation (35) represents the basic input-output rela- 
tions for the amplitude operators in the Fourier space. 
The coefficients therein are determined only by the (com- 
plex) permittivities and thicknesses of the layers of the 
plate and the permittivities of the surrounding media. 
In particular, for s'q n — > (i.e., £ ,« — * 1) Eq. (35) imme- 
diately yields the input-output relations for the special 
case of the plate being surrounded by vacuum. 

The input-output relations in the form of Eq. (35) are 
generally valid, independent of the mechanism of creation 
of the incoming fields and the fields inside the layers of 
the plate. It is worth noting that they take into account 
both propagating waves and evanescent waves. In partic- 
ular in the case when the plate is surrounded by vacuum, 
then the input and output amplitude operators are as- 
sociated with propagating waves (u>/c> k) or evanescent 
waves (u)/c<k) in z-direction. 

The temporal evolution (in the Hciscnberg picture) of 
the amplitude operators is determined by the time de- 
pendence of the basic variables f(r,w) [or equivalently, 
f(z,k, u)], which is governed by the Hamiltonian (2). In 
the special case when the plate is free of active atomic 
sources, then the fields inside the layers of the plate rep- 
resent Langevin noise sources associated with material 
absorption. If absorption is disregarded and the plate is 
surrounded by vacuum that may contain active sources, 
then Eq. (35) solves a scattering problem of the type con- 
sidered in Ref. [23] , namely scattering of fields containing 
evanescent components. In general, absorption cannot 
be disregarded and both atomic sources and Langevin 
noise sources contribute to a field (e.g., in cavity QED), 
which can contain both propagating and evanescent com- 
ponents (e.g., in near-field scanning probe microscopy). 

The input-output relations (35) enable one to calcu- 
late correlation functions of the output field amplitudes 
in terms of those of the input field amplitudes and the 
amplitudes of the fields inside the plate. The simplest 
case is the calculation of the expectation values of the 
field amplitudes. For example, in a typical scattering ar- 
rangement the active light sources are located outside the 
plate, so that the field inside the plate is the absorption- 
assisted random field whose thermal-equilibrium expec- 
tation value vanishes. Application of Eq. (35) thus leads 
to the expectation-value relations 



(E|l(^))=^(Ef»(k,,)) + t: /0 (EW(k,,)), 



(E^ ut (k,a;))=^ /n (E^(k )W ))+< /0 (E^ ( k,a;)) ) 

(44) 

which exactly correspond to standard results in classical 
optics. On the other hand, when the active sources are 
located (in a cavity-like system) inside the plate, we find 



that 

n-l +q 

_ Urii 



3 = 1 



(45) 



e -^(E«(k, W )) + r), e^(E$(k,u;)) 



<E^ t (k,-)H£^e 

. IS an 



(46) 



(E«(k, W ))+r| /0 (E«(k )W )) 



Needless to say that when there are no active sources in- 
side the plate, then (E^2(k, u)) =0 (j = 1, 2, . . . , n - 1) 

is valid and thus (E^j™jj (k, w)) = 0. Clearly, higher-order 
correlation functions of the outgoing amplitude operators 
do not vanish in this case in general. For example, the 
spectral intensity (in the k-space) of the radiation out- 
going from a plane of a plate in thermal equilibrium at 
temperature T, for chosen polarization, is proportional 
to wfoutfcv), where 



(E^t(k,c,)E^(k',-')) 



(47) 



Applying Eq. (35) together with Eqs. (38) and (17) and 
making use of 

= n(uj,T)5^5(uj -w')5(k-k') (48) 
(j ^ 0, n), where 



n{u,T) = 



exp 



\k B T 



- 1 



(49) 



(fee, Boltzmann constant) is the well-known Bose- 
Einstein distribution function, we derive 



«;W t (k, W )=n( W ,T)2 



t q , 



D„ 



\r q iV-?) 
\'j/0\ C Q~ 



,(k,o;)+[r] /oC « + (k, W )+c.c.]} (50) 



(and w^°J ut (k, ui) accordingly), with the coefficients 
(k, uj) (A, A' = ±) being given in Eqs. (A. 15) and 



(43) (A. 16) in the appendix. 



B. Input-output relations in the p-space 

The input-output relations (35) in the k-space can be 
transformed, according to Eq. (15), into the p-space in 
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a straightforward manner. Writing the electric-field op- 
erator E^ ) (r, uj) = E(°) (z, p, uj) at the boundaries of the 
plate as 



we derive 



E<°>(p,a/) = E(°)(«,p,w) 



2=0" 



= E^(p, W )+E^ t (p, W ), (51) 



E( n )(p,w) = EW(«,p,w) 



2=0+ 

J g out V 



EW(p, W )+EW (p,a;), (52) 



E^(p,a;) 



J d 2 p'R 0/n (p,p',uj)--E^(p',uj) + J d 2 p'T n/0 (p,p',uj)--E<£\p',uj) 

n—l . 

f £ J {*$(p, P» ' E + } + *£> (p, p', w) • E« (p' ; W )} . 

u) ■ (p', w) + / d 2 p' T 0/n (p, p', W ) • E| n 0) (p', W ) 



(53) 



J d 2 p' R n/Q (p,p', 

n—l „ 

+ £ J dV {*^(p,^ W )-E«(p', W ) + *«(p,p', W ).E«(p',a;)}, (54) 



where the electric-field operators E^'™^(p, ui) and and 
E2^(p, w) are respectively related to the amplitude op- T /„(p, p',w) 



erators £^° n n) (k, w) and £^ (k, u) (j = 1, 2, . . . , n - 1) as 



= (2^/ d2fc E e^(k)^„e^(k )e -(^'), (60) 

q—p,s 

E L°W) = 7^2 / d2fc E 45(k)ES(k, W ) e ik ", (55) T„ /0 (p,p» 

= (^/ d2fc E eS_(k)tS /0 eW(k)c^), (61) 



(27T) 5 



q=p,s 



q=p,s 



E2'W)= / d 2 ^ E e?2(k)E?i (k, W )e < ^. (57) 



and the integral kernels that relate the outgoing fields to 
the intraplate fields are given by 

*o^(P,P» 



In the above, the tensor-valued reflection and transmis- 
sion integral kernels, respectively, read as 



(2 



*n±(P, P',W) 



y d 2 ^ e^(k)^(k)e^P' 



(62) 



q=p,s 



Ro/n(Pi P',<*>) 

= i/d^E^(kKV^(k)e^ 



= (^/^ E e2W?Leg(k)c*" («S) 



(27T) 

R-n/o(P, P' 
1 

" (2^) 



(58) 



Note that the integral kernels depend only on the dif- 
ference p — p', because of the translational symmetry in 
the x- and y-directions. The input-output relations in 
the form of Eqs. (53), (54) can be considered as a special 
-J d 2 k e^(k)r 9 ^ eg™' ) (k)e ik ( p ~ p '' (59) realization of the general relations given for a dielectric 



q=p,s 



scattering region of arbitrary shape in Ref. [12]. 
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C. Input-output relations in terms of bosonic 
operators in the k-space 

1. General case 

In many applications it may be advantageously to ex- 
press the incoming and the outgoing field operators in 
terms of appropriately chosen bosonic operators. Re- 
calling the basic commutation relations (4) and (5) and 
the orthogonality of polarization unit vectors as given by 
Eqs. (27) and (28), from Eqs. (31) and (32) together with 
Eq. (17) and (18) and Eqs. (36) and (37) we find that the 
input amplitude operators satisfy the commutation rela- 
tions 

kgin l k ' W )' J Vm ^ ' U > 

= cjk n) (k,a;) 8 qg ,6{w - w')<J(k - k'), (64) 

where 

Similarly, the output amplitude operators and the in- 
traplate amplitude operators can be shown to satisfy 
commutation relations of the same type, i.e., 

E(°:)(k lW ),Ejf(k' lW ') 

= 4°^(k,a;)V^-^(k-kO, (67) 

"Eg(k,o;),Ejlt(k', W ')" 

= c%(k,u)5 qq ,8(LO-u')5(k-k>) (68) 

[for the coefficients c q °^l(k, u>) and c 9AA ,(k, w), respec- 
tively, see Eqs. (A.12), (A.13) and Eqs. (A.15), (A.16) 



in the appendix]. Needless to say that input amplitude 
operators, that refer to different sides of the plate, com- 
mute, input amplitude operators and intraplate ampli- 
tude operators commute, and intraplate amplitude oper- 
ators that refer to different layers also commute. Note 
that output amplitude operators that refer to differ- 
ent sides of the plate do not commute in general [see 
Eqs. (A. 9) and (A. 10) in the appendix], which is similar 
as in the one-dimensional situation (cf. [11]). 

From Eqs. (64) and (67) it is seen that bosonic input 
and output operators can be introduced according to 

E£L(k, w) - V / ^t(k^)4° i ' n "L(k, w), (69) 
thus 

« 9 in,out( k ,^),< i„,out( k ^) 

= S qq ,S(uj-uj')S(k-k'). (70) 
Since for A =/= X' in Eq. (68) c^ x , (k, lu) ^ is valid in gen- 
eral, it is useful to introduce intraplate bosonic operators 
according to the superposition 

Eg(k, W ) = ]T r$ A ,(k,a;)ag> (k, W ) (71) 

(j = 1, . . . , n — 1) and choosing the coefficients t^ aa , (k, u) 
[Eqs. (A. 20) and (A. 21) in the appendix] in such a way 
that 

"ag(ic W ),agU(k', w O 

= 5 AA ,V%-w')5(k-k'). (72) 

Substituting Eqs. (69) and (71) into Eq. (35), we may 
express the input-output relations in the k-space in terms 
of bosonic operators, 



a 



<°out(k^) 

out 



f o/n( k ^) *„/o( k >w) 



z (0) 

q in \ 



CV n (k,w) 



(k,w) 



*S/»( k .<") »* 



n/0 



(k,w) 



(k,w) 



n-1 

£ 



%V(k^) 



S-(k,w) 



.0) 

'"</+ 

%0') 



(k,w) 
(k,w) 



(73) 



where the coefficients [modified in comparison with 
Eq. (35)] read 



'0/r. 



(0) 
c • 

(0) '0/n- 
L qout 



t/0 



\ 



(n) 
C ■ 
g m 



f 9 

l 0/r. 



(74) 



f 9 - 

t 0/n ~ 



„(°) 



(n) n/0 
"gout 



f 9 - 
n/0 ~ 



(n) 
r • 

gin q 

>r r «/o ■ 

q out 



(75) 
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and 



10) 

"qO,n A 



(4 ^)" Va E (76) read 



Note that for k = (i.e., perpendicular incidence of light 
on the plate), the input and output relations (73) exactly 
agree with previous results [11] obtained within a one- 
dimensional treatment. 



2. Plate surrounded by vacuum 

Let us turn to the limiting case that the space outside 
the plate - except for possible active atomic sources - 
may be regarded as being vacuum, i.e., 



if lu/c < k, 



y/u> 2 /c 2 - k 2 if w/c > k. 



(77) 



(78) 



For the propagating-field components observed for 

, ..„ , respectively, 



uj/ok, the coefficients cf ; '^ and cf '" 



■*q in 



,(0) 
"gin 



and 



(n) 
u qm 



„(0,n) 
"ij out 



irh 



„(0,n) 
'gin ' 



(79) 



(80) 



as it can be seen from Eqs. (65) and (66) and Eqs. (A. 23) 
and (A. 24) in the appendix. Hence, bosonic input and 
output operators can be defined according to Eq. (69). 
Moreover, the output amplitude operators that refer to 
different sides of the plate commute [Eq. (A. 22) in the 
appendix], which implies that the associated bosonic op- 
erators commute as well. Making use of Eqs. (79) and 
(80), we see that the modified reflection and transmission 
coefficients (74) and (75), respectively, become identical 
with the original ones. For u/c> k, Eq. (73) thus reduces 
to 



a (0) 

u q out 

-(«) 
a q out 



(k,w) 
(k,w) 



*o/„( k > w ) </o( k - w ) 




q in 



(k,w) 
(k,w) 



E 



r 



0$ + (k, w ) <®-(k, W 



^+(k,w) 



n- 



(k,w 




where the transformation matrix connecting the bosonic 
output operators with the bosonic input operators is ex- 



actly the same as that for the corresponding amplitude 
operators in Eq. (35), and the matrix equation 



J 



'0/r. 
t q 

C 0/n 




,J qn+ Vqn- 



J q0+ fqn 
°q0- <Pqn 




(82) 



is valid. In this case, the field operators E' ) (z, k, lo) and Eqs. (29), (30), (36), and (37)] can be represented as 
E(") (z, k, lj) at the boundaries of the multilayer plate [see 



E(°)(0-,k, W ) = -- 1 / 



]T [e^k^k,^) +e(° ) (k)4 o ) ut (k,c) 



q=p,s 



E<»>(0+k, W ) = -^/-f*- 



£ [e^(k)a£ n (k,o;) +eg(k)aiS rt (k, W ) 



(83) 
(84) 



q=p.s 



(/3o=/3n = \A^ 2 /c 2 — fc 2 , Lu/ok). Note that in Eqs. (83) and (84) is nothing said about the location of the active 
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light sources. Clearly, if there are no active light sources 
outside the plate, then E(°'(z,k,w) and E-(")(2,k, lo) 
evolve effectively freely there [cf. Eqs. (33) and (34)]. 

For the evanescent-field components observed for 
lo/c< k, the coefficients cf) n [Eq. (65)] and [Eq. (66)] 
identically vanish, because of (3' — j3' n — 0. Recalling 
Eq. (69), we see that bosonic input operators cannot be 
introduced for the evanescent-field components. Hence, 
it is impossible to extend the validity of Eq. (81) to the 
evanescent-field components. To treat them, one has 
therefore to go back to the generally valid input-output 
relations (35) for the amplitude operators. Clearly, when 
there are no active light sources at any finite distance 
from the plate, then the input field can be regarded as 
being effectively a propagating field [see Eqs. (31) and 
(32)], that is to say, a field that does not contain evanes- 
cent input components. However, the output field may 
contain evanescent output components resulting from the 
intraplate field. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

Applying the Green-tensor formalism of the quantiza- 
tion of the electromagnetic field in the presence of inho- 
mogeneous, causal, linear dielectric bodies, we have de- 
rived the most general three-dimensional quantum input- 
output relations for the electromagnetic-field operators at 
an arbitrary multilayer dielectric plate, and we have given 
the commutation relations needed. Taking into account 
(i) dispersion and absorption of the plate and the sur- 
rounding medium, (ii) possible active light sources inside 
and/or outside the plate, and (iii) both propagating-field 
components and evanescent-field components, the input- 
output relations connect the output field operators at 
the boundaries of the plate with the input field opera- 
tors at the boundaries of the plate and the fields inside 
the layers of the plate. The outgoing fields outside the 
plate can then be obtained from the ones at the bound- 
aries of the plate via quantum Langevin equations, with 
the Langevin noise sources being determined by the ac- 
tive and passive light sources in the respective half-space 
outside the plate. 

We have given the input-output relations in terms of 
both algebraic relations for amplitude operators in the 
k-space and integral relations for field operators in the 
coordinate space. Further, we have shown that in the k- 
space it is possible, provided that the incoming fields con- 
tain effectively propagating components, to rewrite the 
input-output relations in order to obtain them in terms 
of relations between bosonic operators instead of ampli- 
tude operators. Therefore the input-output relations for 
the amplitude operators are more general than those for 
the bosonic operators. 

Finally, we have studied the limiting case of the plate 
being surrounded by vacuum. In this case, the input- 
output relations in terms of bosonic operators only apply 



to the propagating-field components, for which lo j c > k is 
valid. An application to the evanescent-field components 
(iu/c<k) would run into contradictions. To treat them, 
one has therefore to go back to the amplitude operators. 

The derived input-output relations generalize previous 
results obtained from a one-dimensional calculation [11]). 
In the special case of a plate that is embedded in vacuum, 
without active light sources at any finite distance from 
the plate (and without active light inside the plate), they 
reduce to those given in Ref. [21], if evanescent-field com- 
ponents are disregarded. 

In conclusion, the derived input-output relations are 
suited for studying the quantum statistical properties of 
electromagnetic fields in the presence of planar multilay- 
ered structures (such as cavity-like systems or photonic 
crystals), with special emphasis on absorption-assisted 
quantum decoherence, including evanescent-field effects. 
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APPENDIX: COMMUTATION RELATIONS 

Let us begin with the calculation of the commutators 
of the output amplitude operators E^"^ (k, lo). From 
Eq. (35) it follows that 

Ef» t (k, W )=r; n E2(k,,) 

+ ^ /0 E^(k,c,)+F(°)(k,c,) (A.l) 

and 

E^ t (k, W )=^„E;? n (k,a;) 

+ r« /0 Eg(k, W )+FW(k,a;), (A.2) 

where 

n-l 

F(°)(k, W ) = £ [^ + E«(k, W ) + ^_E«(k,u;) 

n— 1 „ 

= ^ e^(k)-V / dz S ^)(0,z,k,c)-j«(z,k,c), 

(A.3) 

n-l 

F(")(k,c) = ]T [^EW( k)W) + ^_E?2(k, W ) 

n— 1 ~ 

= ^ e^(k).V / dzgW\0,z,-k,w)-F\z,-k,w). 

(A.4) 
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Hence, the (relevant) commutators of the output ampli- can be given by 
tude operators that refer to different sides of the plate 



'E(°) ut (k, W ),Efl(k', W ')" 



~ 'o/™ 0/n 



^ i n/0 , n/0 



EW(k, W ),EWt(k>')l + f^(k )W ),^ )t {k',a/ 



= V5(c-c')5(k-k') (4? n r« /n ^;„+cW^ /0 <; ) + [F(°)(k, W ),Ff)t (k>0 



(A.5) 



where Eqs. (64) - (66) have been used. Making use of Eqs. (17) and (18) and recalling the basic commutation relations 
(4) and (5), we derive 



= V^--')5(k-k')^(^) 2 e (° ) ,(k)e^;(k)X: / dzgWH0,z,±u,)£e>!g$>(p,z,t, U ) 

6 j=i J b] 

= V^-cO^k-kO^^'e^^eS^k)/^ / dzg^^zM^'g^^zM 

/0 2 /-oo 2 ~) 

dz<°)(0,z,k^)^4'^ 0) *(0,^k,c,)- / (k ff (°»)(0,«,k ) a;)^^ ) *(0^,k,a;)[. (A.6) 
-oo c JO C 



It is not difficult to calculate the last two integrals the classical Green tensor G(r, r', w), rewritten in terms 

in Eq. (A.6), by using the explicit expression (20) for of g^'fz^'ik, u), 
g^i \z, z',\s., u>). To calculate the sum of integrals 
Sj=o f\j] d^ • • ■, we employ the integral relation (11) for 



£ / d^/'W',^)^ e '> ll gfP\z',z"^) = ±gW\z,z',K<») 

- V. z, k, W ) + </<#'>(*, z', k, W )e, „ % e Zfi , + e xli ^ z, k, w )e a „ . (A.7) 



After lengthy, but straightforward calculations we then 
derive, on using Eqs. (65) and (66), 



Ei u j ut (k, co), (k', J) = 5{u w') S(k k') 



Ei j ut (k,^,E^ u T t (k',c') 



(")t 



E(°) ut (k, W ),Eilt t(k>0 



2 l 0/n 



I A) 



2 n/0 I ' 



(A.9) 



= 0, (A.8) 



12 



7r/i /t^\ 2 
e V c / 



/3„|fc„| 2 Y fc*2 



+ 



r 0/n / 


l,* 2 
l2 k 
i-2 


Po\ko\ 2 \ 


00 ,q* 

|/3 | 2 °/" 







0n ,p 



» 



(«)* 



) - C/o <- ( k ) • <- ( k ) <- (k) • eff -(k) 



(n)* 



(n)* 



e£T(k).eg(k) e^(k).e^(k) U( W - a/) J(k - k') . 



(0)/ 



(A.10) 



Next, let us consider the commutators of the output 
amplitude operators that refer to the same sides of the 
plate. Performing the same steps as before, we now arrive 
at 



Eil t (k, W ),E^ t(k> ')" 

E^dK^^lik'^')] =0, (A.ll) 



Ei°J ut (k, W ), E^k', W ')J = 5(w - u/) <5(k - k') 
irh /w\2 / p 2/3, 



e \cJ 



r s II 
1 1 



\0o\ 2 |/?o| 2 0/ 



(A.12) 



CLt(k,^),E(tt(k>') 



+ 



+ 



'0/r, 



J,* 2 
,2 K 



»| 2 



\0o? 



e^(k)-e^(k) 



7T?l /^\ 2 

e V c) 

k 2 
' N 2 

2 

0/n 



'0/r 



[/3o|A;o| 2 
0o , ft 



l,* 2 



4°j(k).e^(k) 



i- 2 



l2_S_ 



2/3 



IA)| 



|fco| 2 \\0o\ 2 \k o \ 



(0), 



'0/r 



5(w-w')*(k-k / ). 



(A.13) 



The commutators [E^ ut (k,w),E^J lt (k',cj')] are ob- 
tained from Eqs. (A.ll) - (A.13), by making the replace- 



ments fa -> /3n, ^o -> fcnj 4± -> e^ ; , and r^ /n -> r* /Q . 

Finally, it can easily be proved that the intraplate am- 
plitude operators (38) satisfy the commutation relations 
(j = l,...,n-l; A = ±) 

= cg> A ,(k, w) V % - w') 5(k - k'), (A. 14) 



» 



where 



C q±± 



± 



irh /ui\ 
e V c) 



2 P' 



\0A 2 

(e ±2 ^-l) ea(k).eg*(k), 



(A.15) 



Si) 
c q±T 



±i 



.irh /w\ 
eo V c) 



X ( e T2^^ 



2 fL 

\0A 2 

Mi 



l)e^(k).e?>(k). (A.16) 



Using Eqs. (A. 14) - (A.16), we find that the operators 



a?2(k,a;) 



£?2(k>") 



E^(k, W )±E^(k,a;) , (A. 17) 



where 



4±M 



2nL i nHt e -p'/di/2 



C0j V eo 

x {^sinh(/3^-)e«(k).e«*(k) 
±^sin(/3^-)e«(k).e«*(k)} 1/2 



(A.18) 



satisfy bosonic commutation relations. From Eq. (A. 17) 
it then follows that 

£$(k,a;) = £ r« ) v (k, W )ag(k, W ) ) (A.19) 
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where 



and from Eqs. (A. 11) - (A. 13) we find 



7-,+± = ^(k,w)£ 



^ ± =±K«(k,c). 



(A.20) 
(A.21) 



In the limiting case eo.n (when the plate is sur- 
rounded by vacuum) one has to distinguish between 
u)/c>k (propagating-field components) and uj/c<k 
(evanescent- field components). In the first case (uj/c> k) 
we have f3' = /3 = /?„ = f3' n > 0, so that Eqs. (A. 8) - 
(A. 10) reduce to 



E 



(0) 
q out 



(k, W ),Ejl(kV) 



0. 



and Eqs. (A.ll) - (A.13) simplify to 



E 



(0) 



q out 

Ilk /LO\ 2 1 



(k,w),E<?» t (kV) 



g) 2 iv%-"V(k-k'). (A.23) 



E^ ut (k, W ),E^ ut (k>') 



= — (-I' 7T V «(w - w') <5(k - k'). (A.24) 

In the second case (uj/c<k) we have /3' =/3' n = 0. Equa- 
tions (A.8) - (A. 10) then lead to 



E(t t (k,c,),E- ut (k>') 



(")t 



q' out V 



h 2 ^ - J) 5(k - k'), (A.25) 



eo 



l/3o| 



E(°i t (k,c),E^ ut (k',c') 

7Tft /W \ 2 2r 0/n 



p.(°)t 



(c) TffV^-^)^k-k'), (A.26) 



(A. 22) and similarly 



'Eig 1 .(k |W ),E?l(k' )t y); 

= -O 2 ^V^-')^-k0. (A.27) 

eo vc/ |p n | 



Note that the commutation relations (42) in Ref. [21] 
for the bosonic output operators are only valid for prop- 
agating waves, not for evanescent ones. It should also 
be noted that the free-space light modes XJ p .k(z 7 ui) used 
in Ref. [19] are introduced by explicitly requiring that 
u/c> |K|, where K corresponds to k in the present pa- 
per. Hence, these modes solely represent propagating 
waves. 
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